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Abstract 

The double-exchange model for the Mn oxides with orbital degeneracy is 
studied with including on-site Coulomb repulsion, Jahn- Teller (J-T) coupling 
and doping- induced disorder. In the strong interaction limit, it is mapped 
onto a single-band Anderson model, in which all scattering mechanisms can 
be treated on an equal footing. A sharp rise in the mean square fluctuation 
of lattice distortions is found near the Curie temperature T c , in agreement 
with experiments. We show that the spin and J-T disorders lead to a metal- 
insulator transition (MIT) only at low carrier density. The MIT observed in 
samples with 0.2 < x < 0.5 can be explained by further including the disorder 
effect of cation size mismatch. 
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The importance of electron-lattice coupling to the transport properties of the Mn oxides 
Ri-xA x Mn0 3 , which exhibit colossal magnetoresistance (CMR) effect [.lj-f§] in the hole- 
doping range 0.2 < x < 0.5, was pointed out in early theoretical works |5|,§]. Experimental 
measurements |]7|-|T0| have repeatedly given evidences for the existence of Jahn- Teller (J-T) 
effect in the Mn oxides. In those works, the novel ferromagnetic metal (where dp/dT > 
with p as the resistivity) to paramagnetic insulator (dp/dT < 0) transition was understood 
as a crossover from large to self-trapped small polarons caused by the band narrowing 
accompanying with local spin disordering during the transition. Such a description does 
make sense if the carrier density is low enough so that the small polarons are well separated 
in space. In the high density regime, e.g., for the Mn oxides with intermediate doping, there 
should be large overlaps among the polaron wave functions, and the validity of small polaron 
picture becomes questionable. On the other side, recent experiments JIO] showed that large 
and small polarons actually coexist in the insulating paramagnetic phase, an indication that 
even large polarons can also be localized and do not contribute to metallic conductivity. 
These issues suggest that a more accurate criterion for the metal-insulator transition (MIT) 
is urgently needed. 

A valid theory for the MIT at intermediate doping has to be based upon a unified non- 
perturbative treatment of these equally important scattering mechanisms: (i) strong double- 



exchange (DE) interaction between conduction e g electrons and localized spins [11 II}] , (ii) 
J-T coupling and on-site Coulomb repulsion associated with doubly degenerate e g bands, 
and (iii) doping-induced nonmagnetic disorder. Factor (iii) playing an important role in the 
transport properties of the Mn oxides has been well demonstrated in experiments [§,f|]. It 
was proposed there that the size mismatch between R and A atoms leads to bending of local 
Mn-O-Mn bonds, and hence strongly reduces the electron hoping integral. With the basic 
structure of the MnOg octahedra being similar, the compounds Ri_ x A x Mn0 3 with different 
R and A should have about the same J-T coupling strength. Therefore, tunning the cation 
size mismatch has become a practical way to control their transport properties pj. Previous 
theoretical works on the MIT focused on only one or two of these features. For example, in 
Refs. [5] and [6] the lattice effect has been emphasized but the on-site Coulomb interaction 
and doping-induced disorder were neglected. These theories failed to account correctly for 
the doping dependence of the MIT. On the other hand, nonmagnetic disorder effect was 
considered within the single-band DE model in Refs. [15] and [16] with the J-T distortions 
being omitted. The works [15] and [16] overestimated the strength of nonmagnetic disor- 
der necessary for the occurrence of MIT. A complete theory, which is able to describe the 
essential physics of the MIT in the Mn oxides, is still lacking. 

In the present work, the MIT in the Mn oxides will be investigated by incorporating all 
scattering mechanisms outlined above. Employing the path-integral approach, we map our 
Hamiltonian onto a single-band Anderson model in the strong interaction limit. From this 
model, the MIT can be numerically studied using the transfer matrix method without further 
approximation | 17f . The J-T coupling strength is estimated by comparing our calculated 



resistivity with the experimental data of Lai_ x Sr x Mn03. For the first time, the role of the 
J-T distortion on the MIT has been determined very precisely. The phase diagram obtained 
describes consistently the behavior of the MIT observed in samples with CMR as a function 
of doping concentration and cation size mismatch. 

In the Mn oxides, each Mn atom has five outer-shell 3d orbitals, three half-filled t2 9 states 
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giving rise to an S = 3/2 localized spin, and two e g states |+) = \x 2 —y 2 ) and |— ) = \3z 2 — r 2 ) 
forming two-fold degenerate conduction bands. Strong Hund's rule coupling forces electron 
spins on a site to polarize completely, so we can write down a Hamiltonian, in which the 
conduction electrons are effectively spinless, 

H = -J2 fiMiijdj) + E U'n i+ n^ - g Bdjrd,) ■ Q, + ^ £ Q 2 . (1) 

ij i i i 



The first term represents the two-band DE model p3| , |16|| , where dj = (<if + ,<ij_) and 
fij = cos(^/2) cos(0j/2) + sin(^/2) sin^^e - ^ - ^. The second term stands for the on- 
site inter-orbital Coulomb repulsion. The third and fourth terms describe the coupling 
between electrons and two J-T distortion modes Qj = Q x i + Q 2 k, and the harmonic lattice 
deformation energy, respectively ||, where r are the Pauli matrices. The hopping integral 
matrix can be written as [18| 



i ij = (t-r il At)(l + T-n ij ) , (2) 

where t is the bare hopping integral, r^-At is used to model the disorder originated from 
cation size mismatch with At as its amplitude and random numbers between and 
1. Since R and A atoms are randomly distributed, their size mismatch not only decreases 
the average electron band width but also leads to off-diagonal disorder. Here, = n a 
with a = x, y and z for hopping along the x, y and z directions, and n^. = — (v^i + k)/2, 
rij, = (v 7 ^ - k)/2, n z = k. 

Since both the J-T coupling and on-site Coulomb interaction disfavor double occupancy, 
they have the common tendency to split the degenerate e g bands. From Eq. (1), it is easy 
to derive a mean-field Stoner criterion for the occurrence of band splitting at half-filling 
(x = 0) 

(U 1 + 2Ej)D(0) > 1 , (3) 

where Ej = g 2 /k, and -D(O) is the single-band density of states at the band center before the 
band splitting. Equation (3) is also the criterion at which static J-T distortions ((Qi) ^ 0) 
occurs. In the absence of Coulomb repulsion, it needs relatively strong J-T coupling to split 
the e g bands and induce static J-T distortions. In the opposite limit where U' is larger 
than the band width, any small J-T coupling will lead to static J-T distortions. The on- 
site Coulomb repulsion enhances substantially the J-T effect. This is consistent with the 
argument of Varma that the Mn oxides at half-filling (x = 0) are Mott insulators, while the 
J-T effect occurs parasitically . 



According to experiments, the condition in Eq. (3) should be well satisfied in the Mn 
oxides. In order to ensure the system being insulator at x = and to reduce the number of 
involved parameters, we shall project out double occupancy by setting (U' + 2Ej) oo. In 
the path-integral representation, the isospin-charge separation transformation dj = c^yj can 
be introduced, where is the CP 1 field for the isospin and c% is the charge operator [JX9 
The effective Lagrange corresponding to Hamiltonian Eq. (1) is then derived to be 



L ( T ) =J2[ C i C iyi d rYi + 4( d r - V)Ci 



k 

E fiAyliijyMcj - g°i ■ Qi4d + TjQf , (4) 

3 / 
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where Oj(r) = yJ(T)Tyj(r) is a unit vector characterizing the instantaneous isospin orien- 
tation. The third term in Eq. (4) indicates that the isospin kinetic motion is characterized 
by the energy scale K = t(c\cj), which for 0.2 < x < 0.5 is numerically obtained as 
K ~ O.lOt ~ 0.16t at At = 0. Disorder will further decrease this value. On the other hand, 
the J-T coupling given by the fourth term in Eq. (4) is estimated to be —Ej. Therefore, as 
long as Ej > t, the J-T coupling dominates and parallel configuration between Qj and Oj is 
energetically favorable; namely, Oj = Qi/Qi being confined to the x — z plane. A permitted 
parametrization [19 of yj is y] = [cos(0,/2), sin(0j/2)] with — 7r < 0« < vr, where fa is the 



angle between the isospin and the z— axis. 

To study the statistics of the J-T distortions, we can carry out the Gaussian integral of 
the charge variables q(t) in Eq. (4). The saddle point equation of the resulting Lagrange 
yields (Qi) = g(l — x)/k. Then expanding the Lagrange around (Qi) up to quadratic order 
in SQi = Qi — (Qi), we obtain a Gaussian distribution P(e$) = exp(— e 2 /2Aj) /\^2ttAj for 
the dynamical fluctuations of J-T distortions with £j = —gSQi, where 

A 2 j = (k B T)Ej/[l-EjD^)} (5) 

is equal to (ej) = g 2 ((8Qi) 2 )- For fixed Ej, the amplitude Aj/g of lattice fluctuations 
depends on temperature k B T and electron density of states D(ji) at the chemical potential 
H for 5Qi = 0. The denominator of Eq. (5) represents the dynamical aspect of the lattice 
distortions due to electron density fluctuations within the random phase approximation 
(RPA). According to Eq. (4), we obtain an effective single-band Hamiltonian 

H eS = ~ £ tijcjcj + £ i c Ui , (6) 

ij i 

where ty = fij(y\iijYj)i and a constant energy —(1 — x) 2 Ej/2 per site has been omitted. 
With explicit parametrization, £y can be written as 

Uj = fij(t - Tij At) {cos[(fa - fa)/2] + cos[(fa + fa)/2 - la }} , (7) 

with 7 Q = 27r/3, — 27r/3 and for electron hopping in the x, y and z directions. Equation (7) 
is the effective electron hopping integral renormalized by spin disorder, ionic size mismatch 
and orbital disorder. 

To investigate the magnetic and orbital ordering transitions, we replace c\cj in Eq. (6) 
by its average Ko/t. The resulting Hamiltonian for the spins and isospins is studied by 
Monte Carlo simulation on a 20 x 20 x 20 site lattice. We find the spin and isospin ordering 
transition temperatures T c and T*, respectively, to be 

T c ~ 1.3^o , T* ~ 0.5T C . (8) 

If the bare band width W = 12* is taken to be 2eV, then T c is about 250if to A00K for 
0.2 < x < 0.5 in the absence of disorder. The disorder effect may still lower T c by a finite 



amount |]16| . The value of T c evaluated here is comparable to experimental data for the 
Mn oxides. Our finite-size calculation seems to indicate that the isospins are ordered below 
T*, while the result based on a mean- field approximation implies only strong short-range 



ordering [TH] 
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The Hamiltonian Eq. (6) represents an Anderson model with temperature-dependent 
diagonal and off-diagonal disorders. Well below T c , the spins are ordered and the fluctuations 
of the on-site J-T energy e$ are relatively weak, so we expect a metallic phase there. As the 
temperature is increased to above T c , the local spins rapidly become disordered, leading to 
DE off-diagonal disorder. Besides, along with the local spin disordering, the electron band 
narrows and D(/a) increases. From Eq. (5), it follows that the fluctuation amplitude Aj of 
£i also increases near T c . For simplicity, we assume the local spins (6^, <pi) to be completely 
ordered and disordered below and above T c , respectively. The isospins (0j) are disordered 
above T* = 0.5T C according to Eq. (8). If Ej is taken to be Ej = 3.7t (see later discussion), 
by exact diagonalization of Eq. (6) with = on a lattice with 10 x 10 x 10 sites to calculate 
D(fi), it is found that Aj increases by 30% at At = during the magnetic transition. This 
result is consistent with the extensive experimental observations 0HH °f rapidly enhanced 
lattice fluctuations near T c . The increased spin and lattice disorders may possibly lead to 
an Anderson MIT. The picture of Anderson localization is applicable to finite temperatures, 
as long as the electron localization lengths are smaller than the dephasing length due to 
inelastic scattering. This condition could be expected in the Mn oxides in the experimental 
temperature range; otherwise their insulator transport behavior could not be observed. 

Before determining the localization effect due to Eq. (6), we need to estimate the J-T 
coupling Ej(— k 2 /g) by comparing our calculated resistivity with experimental data. Since 
the effect of DE spin disorder and J-T distortions is most prominent in the cleanest systems, 
we consider Lai_ x Sr x Mn03 ||, where the ionic mismatch can be assumed negligible (At ~ 
0). According to the sign of dp/dT, one finds the paramagnetic phase of Lai-^Sr^MnOs 
becomes obviously metallic in the doping range x = 0.3 ~ 0.4 with resistivity of the order of 
10 _3 ncm at T = 400K |J. On the other hand, we can calculate the resistivity from Eq. (6) 
numerically. By using the well-known transfer matrix method developed by Mackinnon and 
Krammer [|17| to calculate the mobility edge, we first determine the phase diagram in the 
n e (= 1 — x) vs A j plane at At = in the paramagnetic phase, as given in Fig. 1(a). Here, 
the disorder coming from both the randomly orienated localized spins and J-T distortions 
has been considered. In the corresponding metallic region, we then apply the Landauer 
formula [^OJ to calculate the conductivity at T > T c as a function of x for different values of 
A j, as plotted in Fig. 1(b), where the experimental data at T = 400-fT of Ref. [9] are shown 
by squares. For 0.2 < x < 0.5, it is found that D(p) changes very slowly with electron 
density n e , so for a given temperature Aj is approximately constant. By comparison of the 
calculated and experimental data in Fig. 2(b), we find Aj~1.5t at T = 400K for the Mn 
oxides, which means Ej ~ 3.71 This value of Ej corresponds to A = 1.1 in the theory of 
Millis et al. H. According to Fig. 1(a), for such strength of J-T coupling, the DE spin 
disorder and J-T distortions lead to an insulator paramagnetic phase only for n e > 0.88 or 
x < 0.12. It is important to notice that since no additional scattering mechanisms besides 
spin disorder and J-T distortions are considered here, the estimated J-T coupling strength 
Ej = 3.7t should be regarded as an upper bound. 

In order to study the MIT for samples with x > 0.12, the effect of the cation size 
mismatch or nonzero At must be included. The electron density of states and mobility edge 
are calculated numerically by finite-size diagonalization and scaling calculation |^TJ. From 



the calculated results, we obtain the phase diagram Fig. 2 in n e vs At plane for the MIT. In 
the limit of large At, the denominator of Eq. (5) may possibly go to zero or negative, giving 
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rise to the RPA instability. Within mean-field approximation, it is easy to see that the RPA 
instability simply means b[ijbn e < 0, an indication of phase separation [22] into hole- rich 
regions with relatively weak lattice distortions and hole-poor regions with strong lattice 
distortions. Such unstable region is also shown in Fig. 2. It is clear that the MIT occurs in 
broad ranges of values of n e and At. With increasing At, localized hole number increases 
and Anderson transition always precedes phase separation. The phase diagram Fig. 2 is in 
agreement with the experimental observation [^,|J that the transport properties of the Mn 
oxides could be drastically altered by tunning the ionic size mismatch. To further determine 
the system properties in the phase-separation regime, one needs to consider anharmonic 
deformation energy of the lattice and take nonuniform states into account. The residue 
resistivity coming from the cation size mismatch in the low-temperature metallic state can 
also be calculated. The critical residual resistivity as shown by the open circles in Fig. 2 for 
occurrence of a MIT is found to be the order of 10~ 4 f2cm or less, being smaller than that 
estimated in the absence of J-T effect ||16|| . The present result seems to be more reasonable, 
and is comparable to data from epitaxial films or single crystal systems ||. 

In summary, employing the path-integral approach and numerical scaling calculations, 
we have shown that the essential physics of the MIT in the Mn oxides with 0.2 < x < 0.5 
can be understood as an Anderson transition driven by DE spin disorder, J-T distortions 
and cation size mismatch. The phase separation predicted in Fig. 2 for strong cation size 
mismatch is also consistent with the inhomogeneous electron states observed recently near 



and above T r in the certain Mn oxides 24 
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FIGURES 



FIG. 1. (a) Phase diagram in n e vs Aj/t plane in the paramagnetic phase where n e = 1 — x is 
the electron density, and (b) the resistivity in the metallic paramagnetic phase as a function of n e 
for several values of Aj/t at T = 400K. 

FIG. 2. Phase diagram for the metal- insulator transition in At/t vs n e plane for Ej = 3.7t. 
In the region labeled "Metal" the system is metallic in both the ferromagnetic and paramag- 
netic phases, while in the regions labeled "M-I(AL)" and "M-I(PS)" the system undergoes a 
metal-insulator transition near T c , where AL or PS indicates Anderson localization or phase sep- 
aration as the origin of the insulator paramagnetic phase. The corresponding critical residual 
resistivity at some points indicated by small circles is shown in unit 10 -4 Slcm. 
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